Universal output directionality of ultrahigh-Q single modes in a deformed microcavity 

with low index of refraction 
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Experimental investigation of the characteristics of quasi-bound states of a quadrupole deformed 
microcavity has revealed five distinct mode groups in cavity emission spectra with cavity quality 
factors different by orders of magnitude and consistently with much different intracavity mode 
distributions but with almost universal far-field emission patterns. These universal directionality of 
high Q modes are explained by a subtle manifestation of unstable manifolds of classical chaos in the 
formation of quasi-bound states. 
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How classical trajectory of a particle is associated with 
its quantum wave function Q is a frequently encountered 
theme in various branches of modern physics, ranging 
from Anderson localization |2j in a disordered conductor 
and its optical version in random media 0, Q to elec- 
tronic wave functions in a high external field in hydrogen 
or in a quantum dot [f| and microwave distribution 
in a Sinai billiard Q. Closely related to this theme, un- 
derstanding mode formation in microscopic cavities is an 
important task from a practical point of view, especially 
in designing efficient laser cavities for photonics appli- 
cations 8]. Particularly, mode formation in deformed 
microcavities of stadium or quadrupole shape has drawn 
much attention since they not only support collimated di- 
rectional output [jj desirable for efficient microlasers but 
also provide a testing ground for classical and quantum 
(wave) chaos [iri|. 

The output directionality observed from deformed mi- 
crocavities has been explained either by ray dynamics 
based on chaos theory || 0] , or by the nature of modes 
obtained from Maxwell's wave equations. Adequacy of 
each approach depends on several factors, but most im- 
portantly on the size of microcavity with respect to the 
wavelength of interest or the size parameter nkr where 
n is the refractive index of the cavity medium, k = 2ir/\ 
the wave vector with A the wavelength and r a represen- 
tative radius of the cavity. 

Most studies on output directionality have been per- 
formed on relatively large cavities with nkr ~ 10 3 or 
so 0, 0, 0] , for which analysis of ray dynamics in 
phase space, corresponding to classical limit of nkr — ► oo, 
has been successful in correctly predicting the output 
emission properties of stadium, quadrupole or hexade- 
capole microcavities Q. 

Although modes of cavity are inevitably involved in 
those experiments with large nkr, measuring actual mode 
spectra and resolving individual modes have not been re- 
ported mainly due to the fact that the free spectral range 



(FSR) of a mode in such large cavities is not much larger 
than the linewidths of individual modes of relatively low 
cavity quality factor Q, associated with the observed out- 
put directionality. 

There have been experiments performed on com- 
parably smaller quardrupole-deformed microcavities 
(QDM's) with nkr ~ 200 |lfj, corresponding to an inter- 
esting regime where wave and particle nature of light may 
coexist. Recently scar modes, i.e. quantum localization 
or a quasi-bound state along an unstable periodic orbit 
(UPO) in ray picture have been observed in deformed 
microcavities with nkr ~ 10 2 - 10 3 [H H H E3 • How- 
ever, spectra of high-Q scar modes with well-defined FSR 
have been only observed in QDM's of low index of refrac- 
tion (n = 1.36) with nkr ~ 200 due to the aforemen- 
tioned technical reasons. 

The aim of the present work is to elucidate the inter- 
play between wave and particle natures of light in mode 
formation in a QDM of properly chosen size parameter 
(nkr ~ 200) with direct comparison between experiment 
and theory with the same size parameter. In particular, 
we report observation of five different well-defined mode 
groups including scar modes, with Q values ranging from 
10 to 10 6 , in both spectrum and output directionality 
in a QDM with a large deformation to enable complete 
classical chaos for the ray undergoing total internal reflec- 
tions. Interestingly, these modes exhibit similar output 
directionality regardless of their different spatial distri- 
butions inside the cavity. This universal directionality of 
high Q modes can be explained by an unexpected man- 
ifestation of ray flow in classical chaos in the formation 
of quasi-bound states in quantum chaos. 

A two-dimensional deformed microcavity used in our 
experiment is made of a liquid jet of ethanol (refractive 
index n=1.361) doped with Rhodamine B dye at a con- 
centration of 10~ 4 M/L. The cross sectional shape of the 
jet is composed of a quadrupole and a much smaller oc- 
tapole component. Nonethless, our cavity can be approx- 
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FIG. 1: Schematic of our experimental setup for measuring 
cavity mode spectrum and single-mode far-field emission pat- 
tern of a QDM. 

imated by a quadrupole, defined as p{4>) = r(l+rj cos 20), 
with a deformation parameter 7/ and a mean radius r of 
about 15 /xm. The details of our new liquid-jet apparatus 
is described elsewhere |l7| . 

As illustrated in Fig. ^ the microjet apparatus is 
mounted in the center of a rotatable stage and the micro- 
jet QDM is illuminated by a cw Ar + pump laser (A p =514 
nm) with its polarization parallel to the QDM column. 
Since the off-resonance pumping efficiency for a QDM de- 
pends on pumping angle |lj| , the pump laser beam is de- 
livered through an optical fiber with its exit end mounted 
on the rotatable stage with a fixed angle of 45° with re- 
spect to the major axis of QDM when output spectrum 
and directionality are measured. The cavity-modified flu- 
orescence (CMF) or lasing light from the QDM is col- 
lected by an objective lens with a full collection angle of 
5 degrees and focused on to an entrance slit of a spec- 
trometer. A polarizer placed in front of the slit selects 
only the polarization component parallel to the QDM 
column. The emission spectrum is then measured for a 
fixed angle 9 of the rotation stage. 

Figure [3 shows the resulting CMF and lasing spectra 
of a QDM with 77=16%. Detection angle 9 was fixed 
at 0=40° for various pumping powers. We can identify 
five distinct mode groups, each of which consists of a se- 
quence of peaks with a well-defined FSR. Modes in the 
same mode group differ from each other by the same 
FSR in frequencies and share similar modal properties 
such as quality factor and spatial distribution. By ex- 
amining their FSR values, we can uniquely label these 
mode groups as I = 1, 2, 3,4, 5 in increasing order of FSR 
in frequency, as shown in the inset of Fig. 2, in analogy 
to radial mode order for a circular cavity. Some mode 
groups undergo laser oscillations as the pumping power 
increases. 

Since the instrumental resolution of spectrometer (0.2 
nm) limits the largest Q value to be directly observed to 
~ 10 3 , we employed two independent methods to mea- 
sure Q values of the observed mode groups. First, in the 
CMF spectra as shown in Fig. the Q value is deter- 
mined from the wavelength at which the visibility of a 
mode drops as we scan from long to short wavelengths 
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FIG. 2: Cavity-modified fluorescence and lasing spectra of 
our QDM with 77=16% for the pump power of 54 mW and 
112 mW. Inset: Five distinct mode groups are identified and 
labeled as I = 1,2,3,4,5 with increasing order of FSR. The 
positions of these mode groups are marked as vertical ticks 
below the spectra. 



since at that wavelength the cavity loss approximately 
equals the known absorption loss of the cavity medium 
[19| . Secondly, the Q value can also be measured from 
the wavelength at which a mode undergoes laser oscilla- 
tion at threshold as seen in Fig. El 16, 20] . Both methods 
consistently give Q = 5 x 10 5 for I = 1. For I = 2, 3, 4, 5 
only the CMF method is used since their Q's are not 
large enough to undergo laser oscillation at the prefixed 
dye concentration. The results are Q = 2 x 10 5 for I = 2, 
Q = 2 x 10 4 for 1=3 and 4 and Q = 1 x 10 4 for I = 5. 

Output emission directionality of each single mode can 
be measured by taking spectra like Fig. [21 for various ro- 
tation angles 9 with respect to the direction of minor axis 
of the QDM and then plotting the peak height of each 
mode as a function of 9, which also equals the far-field 
emission angle in Figs. 3 and 5. By this method, contri- 
bution by background luminescence, mostly due to bulk 
fluorescence of the cavity medium, can be discriminated. 
The resulting single-mode far-field emission pattern is 
shown in Fig. |3] In comparison, previously reported far- 
field patterns for deformed microcavities of nkr ~ 1Q 3 
|l4j ] were measured without any resolution on individual 
modes and thus they correspond to multimode far-field 
emission patterns. 

The most striking feature of the far-field patterns in 
Fig. |21 is that they all look similar, highly directional 
and peaked around 40°-60°. One plausible explanation 
for the observed similarity would be that all five mode 
groups might have very similar intracavity mode distri- 
butions regardless of their completely different Q values. 
However, this is quite unlikely since different Q values by 
many orders of magnitude are usually associated with sig- 



3 




1 ' i ' 1 1 1 ' 

■ ■ 

\ 


1=5 


I 1 1 1 I 1 I 


— -a--- 

r v / V - 

V 1=3 ^jf . 


— affl 

JO 

uniP — mm — ^arcE: 


1=2 


r - m ■ '^-m 

ftri ^ \ 
musfHr:.. ■ * 




1=1 


y - 


-80 -60 -40 -20 





20 40 60 80 



Emission Angle(°) 



FIG. 3: Observed far-field emission pattern of individual sin- 
gle cavity modes. The observed single modes are labeled from 
bottom as I = 1, 2, 3, 4, 5. Baseline for each curve is denoted 
by dotted lines. 

nificantly different mode distributions. In fact, we show 
below theoretically that these mode groups indeed have 
quite different mode distributions inside the cavity. 

Consideration of the ray dynamics in a QDM in clas- 
sical limit (nkr — > oo) often sheds light on mode for- 
mation in a QDM of a finite nkr. The ray dynamics 
can be completely analyzed in terms of a Poincare sur- 
face of section (PSOS) obtained by recording both the 
angle of incidence x an d the azimuthal angle <f> for suc- 
cessive reflections off the cavity boundary. PSOS's for 
a quadrupole with n = 0.16, displayed as backdrops in 
Figs. 01a) and 4(b), show that there exist no stable peri- 
odic orbits above the line of critical angle, sinx c = 1/n. 
None of five mode groups observed in our experiment can 
thus be associated with stable periodic orbits. 

Modes or quasi-bound state solutions of Maxwell's 
equations for an open QDM with 77= 16% and n=1.361 
are calculated by using the boundary element method 
|2ll l2^ |. For direct comparison with experiment, we 
search for modes in the same range of nkr as in the actual 
experiments. In the eigenvalues of the quasi-bound states 
numerically obtained, we also observe five distinct mode 
groups with well-defined mode spacing or FSR, according 
to which modes are labeled as I — 1,2,3,4,5 as before. 
Their Q values are comparable to those of experiment. 

Husimi distributions of some of the modes (/ = 2, 4) 
are shown in Fig. 01 superimposed on the classical PSOS. 
The Husimi distribution is a Gaussian-smoothed version 
of the Wigner function, representing the corresponding 
quantum mechanical probability distribution in phase 
space 0, 0] . We confirm that modes belonging to the 
same mo de g roup show almost the same shape of wave- 
function [25j. Some of the mode groups are localized on 
the UPO's, and thus categorized as scar modes: we can 
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FIG. 4: Calculated Husimi Plots of quasieigen modes with 
(a) nkr — 232.27 and (b) nkr = 226.29, corresponding to 1=2 
and 4 modes, respectively. Distributions below critical line 
(= 1/n) in (a) and (b) are enhanced 100 times in intenisty in 
(c) and (d), respectively. 



associate an octagon with I — 2 mode and a hexagon with 
/ = 4 mode in Fig.^J The five mode groups exhibit quali- 
tatively different spatial distributions, barely overlapping 
in Husimi plot. Complete results on mode calculations 
and analyses on our experiment can be found in Ref . |25l| . 
The far-field emission pattern can be calculated once a 
quasi-bound state is known. Interestingly, all five mode 
groups give rise to very similar far-field patterns as shown 
in Fig. El supporting the experimental observation in Fig. 

m 

So far, it has been widely accepted that the modes 
with distinct shapes exhibit different output directional- 
ity 0, IT3I ITol , so the above experimental and the- 
oretical results contradict to the conventional wisdom. 
This seeming contradiction is, however, resolved by close 
examination on the phase-space mode distribution in the 
Husimi plot. When we enhance the intensity near the 
line of critical angle by 100 times, we can see a few tail- 
like faint structures extending from the regions of wave 
concentration to the region far below the line of criti- 
cal angle. Integrated probability associated with these 
structures below the line of critical angle is negligibly 
small compared to that of the total mode distribution 
for I = 2, 4 in Fig. No matter how small their proba- 
bility is, they are the only direct routes to outside world 
by refraction. Without them, light would have to tunnel 
through a substantial distance in phase space from the 
region of wave concentration to outside across the line 
of critical angle, and thus the probability of such tunnel- 
ing would be extremely low. For the other mode groups 
similar faint structures also exist, providing routes for 
refractive escape of light. 
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FIG. 5: Calculated far-filed distributions associated with 
quasieigenmodes with nfcr=228.87, 232.27, 226.29, 227.65 and 
230.87, which belong to mode groups labeled as I = 1, 2, 3, 4, 5, 
respectively. Curves are smoothed over 5° intervals to match 
the experimental angular resolution. 

Interestingly, the origin of these faint route structures 
in the wave function can be traced back to the classical 
ray chaos. A hint comes from the observation that the 
far-field emission pattern obtained from ray dynamics in 
classical limit for the QDM is similar to those from the 
wave calculation and the experiment. For an open sys- 
tem, long time ray dynamics in chaotic region are pre- 
dominantly determined by the so-called unstable mani- 
folds 0, [yj, which are drawn in Fig. Hfc) and (d) in 
dark lines for our QDM. Since rays escapes from an open 
cavity before reaching completely ergodic limit, ray dy- 
namics is usually restricted in limited phase space and 
thus follows a few dominant unstable manifolds. 

This scenario is, however, valid in the classical limit 
nkr — > oo, and the long time dynamics in this picture 
usually corresponds to several ray circulations, which is 
much shorter than 10 2-3 ray circulations corresponding 
to Q ~ 10 5 - 10 6 observed in our experiment. Neverthe- 
less, the underlying structure of unstable manifold im- 
posed by the geometric shape of the deformed cavity per- 
sists as a weak quantum localization, with 1 /nkr ~ 1/230 
playing a role of an effective Planck constant. Although 
the quantum localization around the unstable manifold 
is not as strong as that for UPO's in the formation of 
scar modes, it leaves a trace as a few faint structures in 
phase-space mode distribution across the line of critical 
angle, resulting in the observed directional output. 

It is noted that for a mode with rather low Q the 
Husimi distribution has significant overlap with the re- 
gion below the critical angle and thus the output direc- 
tionality is mainly determined by the intracavity mode 
distribution as usual |2^|. In addition, for high-Q modes 



for rather small nkr, namely ~ 50, which happens to be 
the size parameter at which many other theoretical stud- 
ies have been performed on the role of unstable manifolds 
in output directionality |l4j |. the faint structure corre- 
sponding to the unstable manifolds has not been observed 
in our numerical studies |25| . 

Our results on universal output directionality can 
greatly simplify design considerations for microcavity de- 
vices; a highly directional output can be achieved regard- 
less of intracavity modal distributions and Q values of 
individual modes in significantly deformed microcavities. 

This work was supported by National Research Lab- 
oratory Grant and by KRF Grant (2005-070-C00058). 
JBS and HWL were supported by a Grant from the Min- 
istry of Science and Technology of Korea. SWK was sup- 
ported by KRF Grant (2004-005-C00044) and by KOSEF 
Grant (R01-2005-000-10678-0). 



* Electronic address: kwan@phya. snu.ac.kr| 
[1] H.-J. Stoeckmann, "Quantum Chaos: an introduction" 

(Cambridge Univ. Press, Cambridge, 1999). 
[2] P. A. Lee and L. V. Ramakarishnan, Rev. Mod. Phys. 

57, 297 (1985). 
[3] D. S. Wiersma, P. Bartolini, A. Lagendijk, and R. Righ- 

ini, Nature 390, 671 (1997). 
[4] A. A. Chabanov, M. Stoytchev, and A. Z. Genack, Nature 

404,850 (2000). 
[5] J. E. Bayfield, G. Casati, I. Guarneri, and D. W. Sokol, 

Phys. Rev. Lett. 63, 364 (1989). 
[6] C. M. Marcus et al, Phys. Rev. Lett. 69, 506 (1992). 
[7] S. Sridhar, Phys. Rev. Lett. 67, 785 (1991). 
[8] K. Vahala, Nature 424,839 (2003). 
[9] J. U. Noeckel and A. D. Stone, Nature 385, 45 (1997). 
[10] L. E. Reichl, The Transition to Chaos in Conservative 
Classical Systems: Quantum Manifestations (Springer- 
Verlag, New York, 1992). 
[11] A. J. Lichtenberg and M. A. Lieberman, Regular and 

Stochastic Motion (Springer- Verlag, New York, 1983). 
[12] C. Gmachl et al, Opt. Lett. 27, 824 (2002). 
[13] T. Harayama et al, Phys. Rev. E 67 015207 (2003). 
[14] H. G. Schwefel et al., J. OPt. Soc. Am. B 21, 923 (2004). 
[15] N. B. Rex et al, Phys. Rev. Lett. 88, 094102 (2002). 
[16] S.-B . Lee et al., Phys. Rev. Lett. 88, 033903 (2002). 
[17] J. Yang et al., Rev. Sci. Instrum. 77, 083103 (2006). 
[18] S.-B. Lee et al, arxiv.org/abs/physics/0603158. 
[19] P. Chylek et al., Opt. Lett. 22, 1723(1991). 
[20] H.-J. Moon, Y-T. Chough, and K. An, Phys. Rev. Lett. 

85, 3161 (2000). 
[21] S. Kagami and I. Fukai, IEEE Trans. Microwave Theory 

and Techniques 32 455 (1984). 
[22] J. Wiersig, J. Opt. A 5 53 (2003). 
[23] Scott Laceyef al, Phys. Rev. Lett. 91, 033902 (2003). 
[24] B. Crespi, G. Perez, and S.-J. Chang, Phys. Rev. E 47, 
986 (1993). 

[25] J.-B. Shim et al, arxiv.org/abs/physics/0603221. 



